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Abstract 

We extend the world model of Kamenshchik et al. to large per- 
turbations by formulating a Zeldovich-like approximation. We sketch 
how this model unifies dark matter with dark energy in a geometric 
setting reminiscent of M-theory. 

After nearly two decades of reign 0], the Einstein-de Sitter dust model 
has been swept aside by observations of high redshift supernovae which 
suggest that the Hubble expansion is accelerating. When combined with the 
Boomerang/Maxima data ^ showing that the location of the first acoustic 
peak in the power spectrum of the microwave background is consistent with 
the inflationary prediction = 1, the evidence for a net equation of state 
of the cosmic fluid lying in the range —1 < w = P/ p < —1/3 is compelling. 
Parametrically, w = Pde/(pde + Pdm) = — gives a ratio of unclustered 
dark energy to clustered dark matter of order 7:3, thereby also resolving the 
longstanding ^dm < 1 puzzle |jl|, Q implied by peculiar velocity fields. The 
theoretical implications of these dicoveries are profound. Simply appending a 
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nonzero cosmological constant A to standard cold dark matter (ACDM, [^) 
invites anthropic arguments as to why both ^dm and Q\ are of order unity 
today. False vacuum models have been proposed , and, in one formulation 
[0], linked to quantum gravity effects and axionic CDM. 

Currently, the most popular alternative is Quintessence which con- 
ventionally entails a real, homogeneous scalar field of mass mg < Hq whose 
potential is arranged such that it silently tracks radiation and matter, grace- 
fully entering a dominant de Sitter phase with an intermediate mixed epoch 
(QCDM) today. Acceptable models can be obtained using pseudo-Goldstone 
bosons , while spintessence |T0| is a twist on this theme involving a complex 



quintessence field. 

If ACDM appears too coincidental, QCDM needs two distinct fields, one 
to describe dark matter, the other dark energy. Economy would suggest 
that dark energy and dark matter should be different manifestations of the 
same entity. Wetterich |Tl| has speculated that dark matter may be cosmon 



(quintessence) lumps while Kasuya has pointed out that spintessence- 
like models are generally unstable to formation of Q-balls which behave as 
pressureless matter. Because they rely upon the coupled nonlinear scalar- 
gravitational fields even into the nonlinear regime, it is difficult to assess 
these proposals in the absence of large-scale numerical simulations. This is in 



marked contrast to I^cdm = 1 where the Zeldovich approximation [|I^ allows 
one to obtain an intuitive picture. Being essentially a variant of ACDM, the 
Barr-Sechel scenario ItII does not suffer this criticism. 



In a recent paper Kamenshchik, Moschella and Pasquier \T^] have studied 
a homogeneous model based on a single fluid obeying the Chaplygin gas 
equation of state 

P = -A/p, (1) 

which has been intensively investigated for its solubility in l+l-dimensional 
space-time, its supersymmetric extension and connection to d-branes [|I^]. 
Implementing this in the relativistic energy conservation equation, the den- 
sity evolves according to 



p{a) = ^A + B/a^, (2) 

where a is the scale factor and B an integration constant. This model 
smoothly interpolates between dust and de Sitter phases without ad hoc 
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assumption. Fabris, Goncalves and de Souza |]T6| have endeavoured to ex- 
amine density perturbations to this model, but, owing to an unfortunate 
choice of the hghtcone gauge, their unperturbed Newtonian equations can- 
not reproduce (0), which is a prerequisite. In any case, the universe is very 
inhomogeneous [|l7l today, and if one takes this idea seriously one needs a 
Zeldovich-like nonperturbative approach adopted to the case P 7^ 0. 

The purpose of this letter is to give such a formulation covariantly and in 
sufficient generality that it can be adopted to any balometric or parametric 
equation of state. Consider the Lagrangian 



(3) 



for a complex scalar field $ = (^/v^?^) exp(— -im^), where m is the mass 
appearing in the potential V . The Lagrangian (^, expressed in terms of 
and 6*, reads 



2n a : ^ 



0,J-K(0V2). 



(4) 



We now apply the Thomas- Fermi approximation ||l8l . In contrast to Kamen- 
shchik et al. [|l^ who assumed spatial homogeneity, we allow for space-time 
variations of 4> on scales larger than m~^. More precisely, we assume that 



(5) 



and hence, we may neglect the second term in the brackets on the right hand 
side of Eq. @ . This yields the Lagrangian 



TF 



(6) 



and the equations of motion for the fields (p and 6 

g'^e^.e,, = v'{<p'/2) , (7) 

(V^0V'^^,.),M = O, (8) 

where V'{x) = dV/dx. Assuming V' > 0, the field 6 may be treated as a 
velocity potential for the fiuid 4-velocity 



(9) 
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that is restricted to the mass shell, i.e., Ufj,U^ = 1. As a consequence, the 
stress-energy tensor T^'^ constructed from the Lagrangian (P) takes the per- 
fect fluid form, with the parametric equation of state 

P = ^V' + V, P = ^V'-V. (10) 



This procedure may be reversed as in the case of Newtonian limit |T9|. Sup- 
pose the equation of state is given, e.g., in parametric form Plip), p{ip)- From 
Eqs. (|TU]) we find 

M^^)=/^. (11) 



P + P 

Using this one obtains p(0^) or ipi^P"^) ciiid hence the potential 

V=l{p-P) (12) 

of the field theory associated with the fluid, the integration constant being 
fixed by the P = limit. For the Chaplygin gas one finds p = (f^ and 



^=oU' + T7 • (13) 




We remark that this potential is very different from that of Kamenshchik et 
al. [0 whose derivation is based on the assumption of homogeneity. Note 
that (0) has the structure of a renormalization group for the fluid-field map. 

For a sensible theory with > w > —1 and a speed of sound satisfying 
< C| = dP/dp < 1, the relativistic limits of these inequalities should 
coincide, which uniquely selects (|T]). It is interesting to note that P = —A/p" 
gives the functional flow 

02(p)=p"(p"+^-A)^ (14) 

which, as p"'^^ — A 0^, resembles that of the periodic gaussian model [^]. 
The critical line cfP' = p describes again the Chaplygin gas and V* = (0^ + 
A/(j)'^)/2 is the corresponding critical potential. Introducing K = cjy^/^fA 
which measures the dimensionless coupling of the 9 field in (|^), the potential 
V* is self-dual under K — > K~^, reminiscent of M-theory string dualities 



. Indeed, the physical content of the full C* shares this duality, as can be 
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seen through ehminating 0^ by its algebraic equation of motion. This yields 



a Born- Infield theory |22| 



C* = -^/AJl-g^-ej^, (15) 



describing space-time as a 3+1 brane in a 4+1 dimensional bulk The 
velocity potential 9 in (^ measures surface excursion in the fifth dimension 
in Eq. (|^. 

Notably, the Chaplygin gas appears in the stabilization of branes in 
Schwarzschild-AdS black hole bulks as a critical theory at the horizon [P^]. 
Also, a free-energy similar to Eq. (|13D occurs in the stringy analysis of black 
holes in three dimensions Replacing Rama's "string bits" with "brane 
cells" , the free energy (pTSD follows on dimensional grounds, corresponding 
to the size of the system. 

Of course the primary question is whether C* yields a reasonable inho- 
mogeneous cosmology. For this purpose we do not eliminate from the 
Lagrangian (0) with (O), writing the field equation (S) in the form 



gp{p + P)U>') =0 (16) 



instead. To solve this equation, it is convenient to use comoving coordinates. 
In the general comoving coordinate system the 4-velocity vector is given by 
= SQ/y/goo- Eq. ( p!6D then becomes 



\p'-A)] =0, (17) 



with the solution 



-^ip' - A) = Bix), (18) 

where B is an arbitrary function independent of x^. For a general metric ^^^j,, 
the proper time is dr = ^/goodx^ and —g/goo = 7 is the determinant of the 
induced 3-metric 

7*i = - 9ij ■ (19) 

9oo 

Noting that for the relevant scales B{x) can be considered to be approxi- 
mately constant, we obtain 

' ^ (20) 
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as the generalization of Eq. (Q), in synchronous coordinates {t = t,x). We 
stress here a pecuharity of negative versus ordinary positive pressure fluids: 
for the former, as the fluid becomes ultrarelativistic, the weak energy condi- 
tion p + P > is saturated so that T^^^-i, = implies = and there are 
no pressure gradients. Thus in both the nonrelativistic and ultrarelativistic 
regimes the synchronous coordinates are comoving. 

Based on this observation, we can take over from the dust case |26] the 
geometric implementation of the Zeldovich approximation as a map from 
Lagrange (q) to Euler (comoving- synchronous, a:;) coordinates inducing the 
3-metric 

7,, = dM'D/. (21) 
Here D is the deformation tensor 

'''' 

and (f the peculiar velocity potential. The function b, which describes the 



evolution of the density contrast, may be calculated in the standard way 
treating the quantity 

h = 2b{t)ip/ (23) 
as a perturbation. From Eqs. (pOD-(p^) one can derive 



with 



where 



p~p(l + 5), P ---(!- 5), (24) 

P 

-p = ^A + B/a\ 5 = ^{l + w), (25) 

w = -4- (26) 

P 

Using these equations and the synchronous metric defined in Eq. (^) with 
we obtain the 0-0 component of the Einstein field equations in the form 

- 3^ + ^/i + ^/i = 47rp[(l + 3w) + (1 - 3w)5]. (27) 

The unperturbed part of this equation 

a An , ^ , , 

- = -—p{l + 3w) (28) 

Qj O 
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Figure 1: Evolution of 6(a)/6(acq) from Ogq = 1.0 x 10 ^ for = 0.7 and 
b'{acq) = 0, for the Chaplygin gas (solid line) and ACDM (dashed line). 

coincides with the second Friedmann equation. The first-order part of (^) 
may be written diferential equation for 6(a), i.e., 

-a%" + a(l - w)b' = (1 + w){l - 3w)b, (29) 
3 

where the prime denotes the derivative with respect to a. The function w{a) 
may be conveniently expressed in terms of the standard parameter Q\ which 
was introduced previously, i.e., 

^(«) = — 7^ 7^ — ^r- (30) 

1 - + ^ao^ 

In Fig. |l| the evolution of b{a) is shown from radiation-matter equality 
with the initial condition 6'(aeq) = 0. Owing to both the high power of a 
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appearing in Eq. (|30| ) and the factor of (1 - 3w) in (^), the difference to the 
standard CDM scenario at a = 1, i.e., today, is neghgable. Note that, while 
the metric perturbation h oc b saturates, the density contrast 5 oc (1 + w)b 
vanishes for large a. In the same figure we give the corresponding result 
for ACDM which is obtained by omitting the factor {l-3w) in Eq. (^) and 
changing to in Eq. (|30|); the contrast 6 in that case is 50% smaller at a 
= 1. Thus for the inhomogeneous Chaplygin gas we can take over the dust 
results bodily up to 2; ~ 0. The picture which emerges is that on caustics, 
where galactic halos and clusters form, we have w ~ 0, i.e., the fluid behaves 
as dark matter. Conversely, in the voids w ~ — ^^a drives the acceleration as 
dark energy. Here the answer to the coincidence question mentioned earlier 
is that acceleration sets in only once the observed cellular structure develops. 

In conclusion, we have shown that the inhomogeneous Chaplygin gas 
offers a simple unified model of dark matter and dark energy. It may be worth 
pointing out that the potential V obtained from the simple equation of state 
(|lD may easily be generalized by adding a more complicated interaction term, 
e.g., a power of greater than 2, and in this way generate a more refined 
equation of state for dark matter, instead of w = 0. This would not alter our 
general picture, as long as this additional interaction term is much smaller 
than \^ when 0^ approaches \/A. In particular, one can add a thermal 
component to (Q) which, although irrelevant to the large-scale structure and 
acceleration, provides a finite phase-space density inferred from galactic cores 
P^ . Ordinary matter can be included using Sundrum's effective field 



theory; the 6 field is derivatively coupled and thus harmless. The fact that 
this unification is achieved in a geometric framework, having roots in the 
"theory of everything" , makes this scenario all the more remarkable. 
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for valuable comments. This research is in part supported by the Founda- 
tion of Fundamental Research (FFR) grant number PHY99-01241 and the 
Research Committee of the University of Cape Town. The work of N.B. is 
supported in part by the Ministry of Science and Technology of the Republic 
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